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TECHNICAL  NOTE  3128 

COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT  FOR 

INTERFERENCE  PRESSURE  FIELD  BETWEEN  WING 

AND  BODY  AT  SUPERSONIC  SPEEDS 

By  William  C.  Pitts,  Jack  N.  Nielsen,  aind 
Maurice  P.  Gionfriddo 

SUMMARY 


Pressure-distribution  data  were  obtained  for  a wing-body  combination 
at  Mach  nvunbers  of  1.48  and  2.00  and  at  Reynolds  numbers  of  0.6,  1.2,  and 
1.5X1C^  to  investigate  the  effects  of  wing -body  interference.  The  model 
was  a single-wedge,  rectangular  wing  mounted  on  a cylindrical  body  with 
an  ogival  nose.  The  body  angle  of  attack  ranged  between  +6°  and  -6°  and 
the  wing -incidence  angle  ranged  from  0°  to  -5.7°.  The  experimental 
pressure -distribution  and  span-loading  results  are  compared  with  the 
linear,  wing -body  interference  theory  of  NACA  TN  2677 . 

For  small  values  of  angle  of  attack  and  wing-incidence  angle  it  was 
found  that  the  experimental  pressure-distribution  results  compared  well 
with  linear  theory,  but  for  larger  angles,  nonlinear  effects  of  angle 
caused  large  differences  from  linear  theory.  The  nonlinear  effects  of 
angle  on  the  wing  were  fairly  well  predicted  by  shock-expansion  theory 
for  the  wing  incidence  case.  In  contrast  with  the  pressure-distribution 
results,  the  lift  loading  was  found  to  be  very  nearly  linearly  dependent 
on  angle.  Reynolds  momber  and  Mach  number  were  found  to  have  only  a small 
effect  on  the  difference  between  experiment  and  linear  theory  except  near 
the  wave  traversing  the  body  from  the  wing-body  Juncture  where  the  effects 
of  both  of  these  parameters  were  large. 


INTRODUCTION 


In  recent  years  much  Interest  has  been  manifested  in  wing-body 
interference.  Some  of  the  theories  that  have  been  developed  for  computing 
the  effects  of  wing-body  interference  on  pressure  distribution  have  been 
compared  by  Phinney  (ref.  l)  and  Lawrence  and  Flax  (ref.  2).  Ferrari 
(ref.  3)  presented  an  iterative  method  based  on  linear  theory.  Morikawa 
(ref.  4)  obtained  sm  approximate  solution  by  solving  a boundary-value 
problem,  and  also  obtained  a closed  solution  by  approximating  the  three- 
dimensional  model  by  a planar  model.  Bolton  Shaw  (ref.  5)  obtained  a 
solution  by  satisfying  boundary  conditions  at  a finite  number  of  points 
rather  than  over  a surface. 
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A method  of  solving  a wide  class  of  wing-body  interference  problems 
is  presented  in  references  6 and  7 together  with  a numerical  application 
to  the  special  case  of  a flat,  rectangular  wing  mounted  at  incidence  on 
a cylindrical  body  at  zero  angle  of  attack.  The  method  determines  an 
interference  potential  that  satisfies  the  boundary  conditions  on  both  the 
wing  and  the  body  when  the  body  is  Introduced  into  the  field  of  the  wing 
alone.  This  interference  potential  is  determined  as  the  sum  of  a number 
of  Fourier  components.  Since  no  assumptions  are  made  beyond  those  of 
linear  theory,  the  exact  linear-theory  solution  could  be  obtained  from 
this  method  by  taking  a sufficient  number  of  terms.  The  mamerical  appli- 
cation shows  that  except  near  the  wing  leading  edge,  the  Fourier  series 
converges  so  rapidly  that,  generally,  only  four  Fourier  components  give 
a close  approximation  to  the  exact  linear -theory  solution. 

Some  experimental  pressure-distribution  data,  such  as  those  of  refer- 
ences 8,  9^  10^  11^  and  12,  are  available  for  comparison  with  theory. 
However,  a complete,  systematic  set  of  data  which  cover  the  effects  of 
angle  of  attack,  wing  incidence,  Reynolds  number,  and  Mach  number  on  the 
pressure  distribution  on  both  the  wing  and  body  of  a wing-body  combination 
is  not  available.  It  la  the  primary  purpose  of  the  present  report  to 
present  such  a set  of  data  and  to  compare  them  with  the  theory  of  refer- 
ences 6 and  7 to  determine  the  accuracy  of  the  theory.  Another  purpose 
is  to  compare  the  calculated  results  of  this  report  with  those  of  other 
theories. 


SYMBOLS 


a body  radius,  in. 

A aspect  ratio  of  wing  formed  by  Joining  exposed  half -wings 

together 

c chord  of  rectangular  wing,  in. 

velocity  amplitude  function  of  n'th  Fourier  component,  in. /sec 
iy  wing-incidence  angle,  radians  except  where  otherwise  designated 

L lift  of  combination  back  to  wing  trailing  edge,  lb 

M free -stream  Mach  number 

n number  of  Fourier  component 

p static  pressure,  Ib/sq  in. 

Pj^  static  pressure  at  any  particular  orifice  of  wing-body  combi- 

nation when  a = ly  = 0,  Ib/sq  in. 
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pressure  coefficient. 
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cylindrical  coordinates: 
(See  fig.  1.) 
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interference  pressure  coefficient  due  to  n'th  Fourier  component 
dynamic  pressure,  Ib/sq  in. 

y = r cos  0,  z = r sin  6 

Reynolds  niomber  based  on  wing-chord  length 
real  part 

semi span  of  wing -body  combination,  in. 
free -stream  velocity,  in. /sec 
characteristic  functions 

Cartesian  coordinates:  x>  axial  coordinate;  y,  lateral  coor- 

dinate; z,  vertical  coordinate,  in. 

(See  fig.  1.) 

body  angle  of  attack,  radians  except  where  otherwise  designated 
wing  angle  of  attack,  deg 

- 1 

dummy  variables  of  integration 

total  wing-alone  perturbation  velocity  potential 

wing-alone  perturbation  velocity  potential  due  to  the  exposed 
right  half  of  the  wing 

wing-alone  perturbation  velocity  potential  due  to  the  exposed 
left  half  of  the  wing 

wing-alone  perturbation  velocity  potential  due  to  the  portion 
of  wing  inside  the  region  occupied  by  body 


Subscripts 


L 

T 


integer  ranging  from  1 to  «> 

lower  surface 
wall  of  wind  tunnel 
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u upper  limit  of  integration 

U upper  surface 

o free-stream  conditions 


EXPERIMENTAL  CONSIDERATIONS 
Apparatus  and  Procedure 


The  investigation  was  made  in  the  Ames  1-  by  3 -foot  supersonic  wind 
tunnel.  This  wind  tunnel  was  equipped  with  a flexible -plate  nozzle  that 
could  be  adjusted  to  give  test-section  Mach  numbers  from  1.2  to  2.2.  The 
pressure  measurements  are  obtained  as  photographic  recordings  of  a 
multiple-tube  manometer  board  using  dibutyl  phthalate  as  the  fluid. 

Since  this  investigation  required  a comparison  of  the  data  for  several 
Mach  numbers  and  Reynolds  numbers  at  the  same  values  of  a and  ly,  it  was 
necessary  to  set  a and  iy  accurately  for  each  measurement  while  \mder 
tunnel  test  conditions.  The  values  of  were  accurately  set  by  means 

of  angle  blocks  in  the  body.  The  angle  of  attack  was  set  by  a special 
image  projection  device.  A mirror  was  inserted  in  the  schlieren  system 
so  that  an  image  of  the  model  was  cast  upon  a screen.  With  the  wind  off, 
the  model  was  set  at  the  desired  value  of  a and  the  inclination  of  the 
model  image  was  marked  on  the  screen.  With  the  tunnel  in  operation  at 
the  desired  pressure,  the  angle  of  attack  of  the  model  was  adjusted  until 
the  inclination  of  its  image  was  parallel  to  the  calibration  line  made  on 
the  screen  with  the  wind  off.  To  check  this  method,  a horizontal  and 
vertical  wire  grid  was  placed  on  the  tunnel  window  and  schlieren  pictures 
were  taken  of  the  model  while  the  tunnel  was  in  operation.  These  pictures 
showed  that  the  image  projection  device  set  a to  within  ±0.07®  of  the 
desired  value.  It  was  especially  necessary  to  set  a accurately  for  the 
small  6ingles  to  avoid  large  percentage  errors  in  the  angle  setting. 

The  model  angle  of  attack  ranged  from  +6®  to  -6°  in  2°  increments 
and  the  wing-incidence  angle  ranged  from  0°  to  -5 >7®  in  1.9®  increments. 

The  test  was  performed  at  the  two  Mach  munbers  1.48  and  2.00  and  at  the 
Reynolds  numbers  of  0.6,  1.2,  ajid  I.5  million,  based  on  the  wing  chord. 

The  model  was  tested  for  all  combinations  of  these  values  of  the  four  para- 
meters investigated. 


Model  and  Support 


The  sting-supported  model,  which  is  diagrammed  in  figure  1,  was  a 
combination  consisting  of  a cylindrical  body  with  an  ogival  nose  and  a 
rectangular,  wedge-shaped  wing.  The  dimensions  of  the  model  are  given  in 
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figure  1.  The  wing  was  made  10  percent  thick  to  minimize  aeroelastic 
effects.  It  was  mounted  in  the  body  by  means  of  a set  of  angle  blocks 
which  enabled  the  flat  wing  surface  containing  the  orifices  to  be  set  at 
0°,  -1.9°>  and  -5 <7°  angles  of  incidence  with  respect  to  the  body 

center  line.  The  pressure  orifices  were  all  located  on  the  upper  surface 
of  the  model.  The  U7  orifices  were  distributed  along  seven  spanwlse  sta- 
tions in  order  to  give  a comparison  with  theory  for  the  wing  and  the  body. 
The  locations  of  the  orifices  are  given  in  table  I. 


Reduction  and  Accuracy  of  Data 


Reduction  of  data. - All  data  are  reduced  to  the  coefficient  form 
(p-p^)  /qQ.  Actually  the  quantity  (p-Pip)/qip  was  measured  and  subsequent 
corrections  were  applied  to  change  the  reference  static  pressure  to  p^^ 
(p^  is  the  static  pressure  at  the  particular  orifice  in  question  when 
a = iy  = 0°)  and  the  reference  dynamic  pressure  to  q^.  Since  p^^ 
Includes  the  effects  of  nose  thickness  and  stream  angle,  using  p^^  as  a 
reference  pressure  minimizes  these  effects  and  essentially  gives  only  the 
pressures  due  to  the  angle  settings  of  the  model.  The  dynamic  pressure 
was  adjusted  from  q^  to  q^  on  the  basis  of  a previous  pressure  survey 
of  the  tunnel.  This  later  adjustment  was  negligible  for  M = 1.48  and 
amounted  to  less  than  a 3“percent  correction  for  M = 2.00.  For  the  pur- 
pose of  comparison  with  theory  the  pressure  coefficient  (p-Pj^)/qp  is 
reduced  to  the  parameters  pP/a  for  iy  = 0°  and  pP/iy  for  a = 0°. 

Accuracy  of  data.-  There  are  two  types  of  errors  that  enter  into 
an  experimental  Investigation:  systematic  errors  and  random  errors.  In 

this  paper  accuracy  will  be  taken  as  the  ability  of  the  experiment  to 
give  the  true  values  without  nose  effect  or  stream  angle  and,  hence,  is 
a measure  of  the  systematic  errors.  Precision  will  be  taken  as  the  abil- 
ity to  repeat  the  data  and,  hence,  is  a measure  of  the  random  errors  in 
the  experiment. 

Several  factors  contributed  random  errors.  The  major  factor  was 
the  error  in  the  angle -of -attack  setting.  The  uncertainty  in  each  angle 
setting  was  t0.07°,  but  each  measurement  was  dependent  upon  two  angle 
settings:  the  setting  for  the  condition  represented  and  the  setting  to 

determine  the  zero  correction.  This  leads  to  a net  uncertainty  of  0.1° 
which  would  account  for  a 5-percent  error  for  angles  of  ±2°.  Most  of 
the  remainder  of  the  uncertainty  in  the  data  is  due  to  the  fact  that  the 
reference  wall  static  pressure  in  the  tunnel  changed  slightly  from  run  to 
run  while  the  total  pressure  remained  constant.  Although  the  magnitude 
of  this  pressure  change  was  quite  small,  it  was  large  enough  compared  to 
the  small  pressure  differences  for  the  2°  angle  settings  to  cause  as  much 
as  a 3-percent  error.  In  addition  to  these  factors,  between  1-percent 
and  2 -percent  xincertainty  was  observed  in  reading  the  data  from  the  mano- 
meter board  pictures. 
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To  determine  experimentally  the  precision  of  the  data,  a large  number 
of  repeat  measurements  were  taken  and  compared.  It  was  found  that  for  a 
or  iy  = ±2°,  two  independent  determinations  of  pp/a  or  pP/iy  differed 
from  each  other  by  ±7  percent  on  the  average . For  a or  1^  = ± and  a 
or  iu  = ±6°,  the  experimentally  determined  precision  of  pP/a  and  PP/ly 
is  ±4  percent  and  ±2  percent,  respectively.  The  precision  in  pP/a 
increases  with  the  magnitude  of  the  angle  because  a large  part  of  the  ran- 
dom error  is  due  to  the  angle  setting.  The  known  major  experimental 
errors  are  due  to  stream  angle  and  body  nose  effects.  The  effect  of 
these  factors  was  not  determined,  but,  as  described  in  the  reduction  of 
data  section,  corrections  were  applied  to  minimize  their  effect,  assuming 
the  effects  did  not  vary  appreciably  with  angle -of -attack  settings.  This 
assxjmptlon  should  be  good  for  the  body-thickness  effect.  However,  it  is 
not  necessarily  a good  assumption  for  the  stream-angle  effect  since  the 
stream  angle  varies  slightly  at  different  positions  in  the  tunnel  and  the 
model  moves  approximately  6 inches  in  a vertical  direction  between 
a = +6°  and  a = -6°.  Since  the  stream -angle  correction  that  was  used 
was  obtained  for  the  a = 0°  position  in  the  tunnel,  data  obtained  at 
a = 0°  should  have  no  appreciable  error  due  to  stream  angle.  For  other 
values  of  a,  some  error  due  to  stream  angle  is  possible.^ 


For  the  purposes  of  this  paper,  the  important  question  is,  "How  well 
does  theory  predict  the  experimental  data?"  Direct  comparisons  between 
linear  theory  and  experiment  will  be  made  only  for  a = ±2°  and  ly  = -1.9° 
data.  In  figure  2,  experimental  pressure  distributions  in  the  wing -body 
juncture  obtained  from  two  independent  measurements  with  iy  = -1.9°  and 
a = 0°  are  shown  together  with  a faired  curve  of  their  average  values. 

The  ±7 -percent  limit  of  precision  about  the  average  value  is  represented 
by  the  dotted  lines.  The  figure  shows  that  the  theoretical  value  generally 
lies  between  these  dotted  lines  and  therefore  the  theory  predicts  the 
experimental  values  within  the  precision  of  the  data  in  this  example. 


THEORETICAL  CONSIDERATIONS 


The  theory  of  this  report  is  a direct  application  of  the  general 
linear  theory  of  references  6 and  7 to  two  special  conditions  of  a rec- 
tangular wing  and  body  combination.  The  mathematical  model  consists  of 
a rectangular,  flat-plate  wing  mounted  on  an  infinitely  long  circular- 
cylinder  body.  The  aspect  ratio  of  the  wing  is  large  enough  so  that  the 
Mach  lines  from  the  leading  edge  of  the  wing  tips  do  not  Intersect  the 
wing-body  juncture.  The  two  special  conditions  treated  are  the  angle -of - 
attack  case  (variable  angle  of  attack  with  iy  = 0°)  and  the  wing  incidence 
case  (variable  wing-incidence  angle  with  a = 0°).  A detailed  example 
calculation  of  the  wing -incidence  case  is  given  in  reference  7>  so  it  will 

^A  stream-angle  and  pressure  survey  of  the  wind  tunnel  indicated  that 
stream-angle  variation  caused  the  magnitude  of  the  experimental  values 
of  pp/a  to  be  k percent  high  on  the  average. 
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not  be  discussed  in  detail  here.  The  calculations  involved  in  the  angle- 
of -attack  case  will  be  outlined  in  the  following  section.  The  differences 
from  linear  theory  due  to  finite  values  of  a and  iy  are  predicted  by 
shock-expansion  theory. 

In  the  analysis  of  references  6 and  7 and  in  the  calculations  for 
this  report  the  body  radius  is  taken  as  unity,  and  is  taken  as  2 so 

that  p = vNr  - 1 = 1.  Any  formula  can  be  generalized  to  any  body  radius 
by  dividing  all  length  symbols  by  a and  can  be  generalized  to  any  Mach 
number  by  dividing  all  streamwlse  lengths  by  p,  by  multiplying  all  pres- 
sure coefficients  by  p,  and  leaving  all  potentials  and  span  loadings 
unaltered . 


Sinnmary  of  Method  for  Angle-of -Attack  Case 


As  in  all  applications  of  the  theory  of  references  6 and  7^  the  velo- 
city potential  of  the  combination  is  considered  to  be  the  sum  of  the  wing- 
alone  potential  and  an  interference  potential.  Since  the  wing-alone  poten- 
tial can  be  determined  by  existing  methods,  the  essential  problem  is  to 
determine  the  interference  potential.  This  is  done  by  developing  a poten- 
tial that  cancels  the  flow,  due  to  the  wing  alone,  across  the  body  surface 
as  well  as  satisfying  the  two  other  boimdary  conditions,  (l)  that  it  does 
not  distort  the  shape  of  the  wing  and  (2)  that  it  is  zero  ahead  of  the  wing 
leading  edge. 

The  first  step  in  calcvilating  the  pressure  coefficients  is  to  deter- 
mine the  wing-alone  potential.  From  this,  the  normal  velocity,  and  hence 
the  boundary  condition  on  the  body  surface,  is  obtained.  The  expression 
for  the  normal  velocity  is  expanded  in  a Fourier  series  of  the  form 

«s 

^2^  ^an(*)  where  TgnCx)  is  called  the  velocity  amplitude  fvuic- 

tion.  After  the  velocity  amplitude  functions  are  determined,  the  inter- 
ference pressure  coefficient  at  any  point  in  the  flow  field  can  be 
directly  found  in  series  form  with  the  aid  of  the  universal  Wan(x,r) 
functions  developed  and  tabulated  in  references  6 and  7 Tor  r = 1.  By 
adding  the  interference  pressures  to  the  wing-alone  pressures,  the  pres- 
sure coefficients  for  the  wing-body  combination  are  obtained.  The 
details  of  the  calculation  for  the  angle -of -attack  case  are  presented  in 
Appendixes  A,  B,  and  C.  It  is  shown  that  the  interference  pressure 
coefficient  for  r * a is 
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Wing-Incidence  Case 


The  theoretical  values  for  the  vlng-incidence  case  used  in  this  report 
are  taken  directly  from  reference  ?•  The  only  exception  is  the  pressure 
distribution  on  the  0 = n/4  meridian  on  the  body  which  is  not  presented 
in  reference  ?•  However,  it  was  calculated  in  the  same  manner  as  was  used 
for  the  top  meridian  of  the  body  so  no  further  discussion  is  necessary 
here. 


RESULTS  AMD  DISCUSSION 


A complete  set  of  data  in  the  form  of  P for  the  Reynolds  niimbers 
0.6,  1.2,  and  I.5XIO®  at  M = 1.48  and  for  R = 1.5X10®  at  M = 2.00  is 
presented  in  table  II.  These  values  of  P are,  for  the  most  part,  aver- 
ages of  two  readings. 


General  Physical  Principles 


Before  discussing  the  results  of  the  investigation  in  detail,  it  is 
well  to  give  first  a general  physical  description  of  the  effects  to  be 
expected.  Figures  3 and  4 show  qualitatively  the  pressure  distributions 
to  be  expected  on  a rectangular  wing  and  body  combination  for  the  angle -of - 
attack  case  and  the  wing-incidence  case,  respectively.  The  chordwise 
variations  of  the  coefficient,  pp/a  or  pP/iy,  are  shown  for  five  stations 
by  the  shaded  areas The  variation  in  the  pressure  coefficient  is  some- 
what exaggerated  for  emphasis.  These  figures  show  that  Mach  cones  emanat- 
ing from  the  wing-body  Juncture  determine  the  points  at  which  the  various 
effects  of  wing-body  Interference  are  felt.  On  the  cylindrical  body  the 
pressure  coefficient  is  zero  in  front  of  the  Mach  helix  originating  at  the 
leading  edge  of  the  wing -body  Juncture.  However,  as  shown  by  the  two  sta- 
tions on  the  body,  the  pressure  rises  abruptly  behind  this  Mach  helix, 
point  1,  in  both  figures . The  Mach  helices  from  the  two  wing  panels  cross 
the  6 = n/2  station  simultaneously  so  that  there  is  only  one  large 
increase  in  the  magnitude  of  the  pressure  coefficient.  These  Mach  helices 
cross  the  6 = 3«/4  station  at  two  different  points  so  that  beyond  point  1 
there  is  a secondary  increase  in  the  pressure  coefficients  at  point  2. 

These  Mach  helices  continue  to  curl  around  the  body  until  they  strike  the 
wing  panel  at  points  3,  where  part  of  the  pressure  disturbance  continues 
along  the  wing  and  part  of  it  is  reflected  along  another  Mach  helix  on  the 

^At  y/a  =3  in  figure  3 the  shape  of  the  curve  behind  the  Mach  wave  has 
been  inferred  from  the  wing-incidence-case  calculations. 

The  pressure  distribution  shown  for  the  0 = 3jt/4  station  on  the 
body  is  identical  to  the  pressure  distribution  for  the  0 = jt/4  station 
due  to  the  symmetry  of  the  model. 
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body,  causing  a further  increase  in  the  magnitude  of  the  pressure  coef- 
ficients at  points  4.  Another  pressure  disturbance  originates  at  the 
trailing  edge  of  the  wing-body  Juncture  that  causes  the  decrease  in  the 
magnitude  of  the  pressure  coefficients  noted  at  points  5 of  the  two 
figures . 

On  the  wing  of  the  combination  the  pressure  coefficient  is  the  same 
as  that  for  a wing  alone  in  front  of  the  Mach  wave  from  the  wing-body 
juncture,  except  that  when  the  body  is  at  an  angle  of  attack  the  body 
upwash  effectively  twists  the  wing  in  a manner  such  that  oy  = a(l+a^/y^) . 
Figure  3 shows  this  effect  of  body  upwash  along  the  leading  edge  of  the 
wing  where  the  pressure  coefficient  decreases  as  y/a  increases  because 
of  the  effective  twist  of  the  wing.  The  importance  of  body  upwash  can  be 
seen  by  comparing  the  pressure  distribution  along  the  leading  edge  in 
figure  3 with  that  in  figure  4.  The  pressure  coefficient  at  the  wing-body 
juncture  in  figure  3 is  twice  that  in  figure  4 where  there  is  no  body 
upwash.  The  pressure  coefficient  at  any  given  spanwise  station  remains 
nearly  constant  between  the  wing  leading  edge  and  the  Mach  wave  from  the 
wing -body  juncture-  Behind  the  Mach  wave,  interference  from  the  wing -body 
juncture  causes  the  pressure  coefficient  to  decrease  in  magnitude  as  shown 
in  the  two  figures . 


Effects  of  Angle  of  Attack 


Comparisons  between  theory  and  experiment  for  the  angle-of -attack  case 
are  made  in  figure  5 for  data  at  a Reynolds  nxamber  of  1.5x10®  and  Mach  num- 
bers of  1.48  and  2.00  with  iy  = 0°  and  a = ±2°  and  ±6°. 

Pressure  distribution  in  J\xncture  of  wing-body  combination*-  A com- 
pari son  between  linear  theory  and  experiment  for  the  pressure  dl s tribution 
in  the  wing-body  junct\a^  is  made  in  figures  5(a)  and  5(b)  for  both  Mach 
numbers.  The  sketches  show  the  pertinent  Mach  lines  and  the  spanwise 
location  of  the  orifices.®  The  experimental  data  points  from  the  wing 
s\irface  on  which  a compression  occurs  (negative  angle  of  attack)  are 
represented  by  flagged  symbols,  and  the  data  points  from  the  surface  on 
which  an  expansion  occurs  (positive  angle  of  attack)  are  represented  by 
unflagged  symbols.  The  figures  show  that  the  theory  predicts  the  magni- 
tude of  pP/a  about  5 percent  below  the  average  of  the  a = ±2°  experi- 
menteO.  values  at  M * 1.48  and  about  15  percent  below  experimental  values 
at  M = 2.00.  The  chordwise  variation  is  well  predicted  by  the  theory. 

Linear  theory  predicts  that  the  parameter  3P/a  is  independent  of 
angle  of  attack.  Actually  it  is  not,  and  the  nonlinear  effects  of  angle 
of  attack  cause  a spread  in  the  data.  The  magnitude  of  the  spread  to  be 

®The  location  of  these  Mach  lines  is  only  qualitative  because  the  calcu- 
latlons  were  made  using  shock-expansion  theory,  with  the  assvnnptlon  that 
there  was  no  local  Mach  number  variation  behind  the  leading  edge  of  the 
wing.  To  simplify  the  sketches,  the  Mach  helices  on  the  body  are  repre- 
sented as  straight  lines.  
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expected  between  a = +6°  and  a = -6°  was  calculated  at  the  leading  edge 
of  the  wing  by  shock-expansion  theory.  As  in  the  calculation  of  first- 
order  effects,  the  local  angle  of  attack  was  determined  using  Beskin 
upwash  theory.  For  M = 1.48  body  upwash  caused  the  shock  wave  to  be 
detached  from  the  wing  in  the  wing-body  juncture  so  that  no  calaHation 
of  the  spread  co\ild  be  made  there.  For  M = 2.00  it  was  found  that  near 
the  wing-body  juncture  the  predicted  spread  was  about  twice  the  experi- 
mental spread,  whereas  for  y/a  greater  than  about  1.5  the  experimental 
spread  was  fairly  well  predicted.  This  difference  between  shock-expansion 
theory  and  the  experimental  data  in  the  wing-body  juncture  is  probably 
due  to  the  combination  of  two  things.  First,  near  the  wing -body  juncture 
the  body  upwash  is  modified  by  viscous  effects.  Second,  the  theoretical 
spread  was  calculated  at  the  leading  edge  of  the  wing  and  this  value  was 
assumed  to  apply  rearward  to  the  first  orifice.  This  assumption  is  prob- 
ably good  beyond  y/a  =1.5  where  the  chordwise  changes  in  pressure  are 
small  back  to  the  first  orifice,  but,  in  the  juncture,  the  changes  in  the 
chordwise  direction  are  large  near  the  wing  leading  edge  so  that  this 
assumption  is  probably  invalid. 

Another  phenomenon  not  predicted  by  linear  theory  is  shown  by  figure 
5(a) . The  linear  theory  predicts  that  the  Mach  helix  from  the  opposite 
wing  panel  (see  sketch)  should  intersect  the  wing -body  juncture  at  point 
1,  causing  an  Increase  in  the  magnitude  of  pp/a.  This  effect  is  observed 
experimentally  for  negative  values  of  a but  in  front  of  point  1 rather 
than  exactly  at  point  1.  The  reason  is  that  for  negative  values  of  a 
a compression  occurs  on  the  orificed  surface  reducing  the  local  Mach  num- 
ber from  the  free-stream  Mach  number,  thus  increasing  the  Mach  angle  and 
causing  the  Mach  helix  to  shift  forward.  The  result  is  the  spread  of  the 
data  shown  in  figure  5(a)  near  point  1.  This  effect  is  not  shown  by 
figure  5(^)  because  the  wing  chord  is  effectively  shorter  for  M = 2.00 
so  that  the  orifices  do  not  extend  to  the  Mach  helix  as.  shown  by  the 
sketch. 

Figures  5(a)  and  5(^)  show  that  Mach  number  has  no  effect  upon  the 
magnitude  of  the  higher-order  spread  due  to  angle  of  attack  or  upon  the 
chordwise  variation  of  pp/a,  but  on  the  average  the  magnitude  of  pP/a 
is  about  10  percent  higher  for  M = 2.00  than  for  M = 1.48. 

Pressure  distribution  on  top  meridian  of  body  of  wing-body  combina- 
tion.- A comparison  between  the  linear  theory  and  experiment  for  the  pres- 
sure distribution  on  the  top  meridian  of  the  body  is  made  in  figures  5(c) 
and  5(d) . These  figures  show  that  theory  and  experiment  are  in  good 
accord  for  a = ±2°,  particularly  at  M = 1.48.  However,  nonlinear 
effects  due  to  a cause  a large  spread  between  the  data  for  a = +6°  and 
a = -6°.  All  the  effects  predicted  to  occur  on  the  body  in  the  section 
of  the  report  "General  Physical  Principles"  are  observed  experimentally, 
but  not  exactly  at  the  points  predicted  because  of  nonlinear  effects.  The 
pressure  rise  predicted  at  point  1 of  figures  5(c)  and  5(d)  occurs  pre- 
maturely and  is  less  abrupt  than  expected  for  all  angles  of  attack  because 
of  the  boundary  layer  on  the  body.  The  variation  in  local  Mach  niomber 
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causes  the  Mach  helices  to  shift  forward  for  the  negative  angles  of  attack 
as  discussed  in  the  section  treating  the  wing -body  Juncture.  The  increase 
in  the  magnitude  of  pp/a  expected  at  point  2,  x/pa  = 3n/2,  actually 
occurs  at  about  x/pa  = 4 for  a = -2°.  The  decrease  in  magnitude  of  pP/a 
that  is  expected  at  point  3 actually  occurs  at  about  x/pa  = 4.0  for 
a = -6°.  For  the  positive  angles  of  attack  the  Mach  helices  are  shifted 
rearward  so  that  these  effects  are  not  observed  experimentally  in  the 
range  of  x/pa  measured. 

Figures  5(c)  and  5(d)  show  that,  in  general,  the  M = 1.48  data  are 
predicted  better  by  the  theory  than  are  the  M = 2.00  data.  For  M =2.00 
there  is  an  unexpectedly  large  pressure  coefficient  in  front  of  point  1 
for  negative  angles  of  attack."*  For  a = -2°  and  M = 2.00,  pP/a  dips 
slightly  near  point  1 and  then  rises  and  overshoots  the  a = -6°  data. 

This  effect  is  due  to  the  boundary-layer  condition  on  the  body  and  will 
be  discussed  in  detail  in  the  section  dealing  with  Reynolds  number  effect. 

Pressure  distribution  on  Q = 45°  meridian  of  body.-  A comparison 
between  the  linear  theory  and  experiment  for  the  pressure  distribution 
on  the  B = 45°  meridian  of  the  body  is  made  in  figures  5(c)  and  5(i’)* 
Essentially  the  same  effects  are  shown  on  this  meridian  as  on  the  top 
meridian. 

Just  as  for  the  top  meridian  of  the  body  the  experiment  is,  in  gen- 
eral, better  predicted  by  the  theory  for  M = 1.48  than  for  M = 2.00, 
and  the  same  boundary- layer  effects  are  evident  near  point  1 for  M = 2 .00. 

Pressure  distribution  on  wing  of  wing-body  combination.-  Experimental 
chordwise  pressure  distributions  on  the  wing  are  shown  in  figures  5(g)  to 
5(n)  for  the  four  spanwise  orifice  stations  y/a  = 1.25,  1»92,  2.58^  and 
3.92.  In  front  of  the  Mach  cone  from  the  wing -body  Juncture  no  interfer- 
ence is  felt  from  the  wing-body  Juncture  so  that  the  theoretical  pressure 
distribution  for  a wing  alone  in  the  body  upwash  field  is  used  in  this 
region  (see  Appendix  C).  The  theoretical  results  behind  the  Mach  cone 
from  the  wing-body  Juncture  (the  region  beyond  point  2 in  the  figures)  are 
not  available  because  the  W2^(x,r)  functions  for  r > a have  not  been 
calculated.  Figures  5(g)  "to  5(n)  show  that,  in  general,  the  wing-alone 
theory  predicts  magnitudes  of  pP/a  about  5 percent  below  the  measure- 
ments for  a = ±2°  for  M = 1.48  and  about  12  percent  below  the  measure- 
ments for  M = 2.00.  The  spread  in  the  data  between  a = +6  and  a = -6 
is  fairly  well  predicted  by  shock-expansion  theory  for  y/a  greater  than 
about  1.5  (figs.  5(1)  to  5(n)).  At  y/a  = 1.25  the  predicted  spread  is 
too  large.  Just  as  was  the  case  for  the  wing-body  Juncture. 

Some  of  the  Interference  effects  discussed  in  the  section  of  the 
report  entitled  "General  Physical  Principles"  are  illustrated  in  figures 
5(g)  to  5(n) . The  interference  effect  from  the  opposite  wing  panel  is 
observed  in  figure  5(g)  where.  Just  in  front  of  point  1,  the  same  spread 
^The  possibility  that  this  large  pressure  coefficient  was  due  to  body 
crossflow  was  considered  by  including  the  second-order  terms  in 
Bernoulli's  equation  in  the  calculation  of  the  body-alone  pressure  coef- 
flclents,  but  this  predicted  only  a small  part  of  the  observed  values. 
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in  the  data  occurs  as  in  the  wing -body  Juncture.  According  to  linear 
theory  the  disturbance  originating  at  the  nearer  wing -body  Juncture  should 
be  felt  at  point  2 of  figures  5(i)  'to  5(m)>  and  the  magnitude  of  pP/a 
should  begin  to  decrease  from  the  wing-alone  value  there.  These  figures 
show  that  the  magnitude  of  pP/a  does  decrease  in  the  neighborhood  of 
point  2.  They  also  show  that,  in  general,  the  a = +6°  and  the  a = -6° 
data  converge  in  the  neighborhood  of  point  2-  This  convergence  is  due  to 
a variation  in  the  local  Mach  number  with  a.  This  is  shown  by  the  sketch 
in  figure  5(j)  where  the  disturbance  from  the  wing-body  Juncture  is  first 
felt  at  point  3 for  a = -6°,  whereas  it  is  first  felt  at  point  4 for 
a = +6^.  Since  the  magnitude  of  pP/a  begins  to  decrease  as  soon  as  this 
disturbance  is  felt,  the  magnitude  of  pP/a  begins  to  decrease  at  a 
smaller  value  of  x/pa  for  a = -6°  than  for  a = +6°,  thus  causing  the 
convergence  observed.  The  sketches  in  figures  5(h)  and  5(®)  show  that  the 
disturbance  from  the  wing  tip  should  also  cause  the  a = +6°  and  a = -6° 
data  to  converge  beyond  point  6 in  these  figures.  The  figures  show  that 
the  data  not  only  converge  but  actually  cross  over  and  reverse  order  Just 
beyond  point  6. 

The  only  significant  effect  of  Mach  number  shown  by  figures  5(g)  "to 
5(n)  is  the  approximately  10-percent-larger  values  of  pP/a  for  M = 2.00 
than  for  M = 1.48.  Nearly  40  percent  of  this  difference  may  be  due  to 
differences  in  stream  angle  in  the  wind  tunnel  for  the  two  Mach  numbers. 

Span  loading  distribution.-  Span'  loading  is  defined  for  both  the  body 
and  the  wing  as  the  integral 


The  experimental  results  for  the  span  loading  distribution  on  the  wing  and 
body  of  the  combination  are  presented  in  figure  6 and  a comparison  is  made 
with  theory  on  the  body.  As  has  been  mentioned,  the  theoretical  results 
on  the  wing  in  the  interference  region  between  the  body  and  point  1 were 
not  available  for  the  a case  at  the  time  this  paper  was  written.  The 
theoretical  pressure-distribution  curves  shown  beyond  point  1 are  those 
of  a wing  alone  with  no  body  upwash  present  since  outboard  of  this  point 
its  effect  is  small.  The  decrease  in  span  loading  beyond  point  2 due  to 
wing-tip  interference  is  calculated  by  Busemann  theory,  reference  13. 
Figure  6 shows  that  theory,  on  the  average,  predicts  the  experimental  span 
loading  on  the  body  within  10  percent.  Of  particular  interest  is  the  fact 
that,  in  general,  the  higher-order  differences  due  to  a that  were  so 
large  for  the  pressure-distribution  results  are  negligible  for  the  span 
loading  distribution.  The  only  exception  is  on  the  top  of  the  body, 
y/a  = 0,  at  M = 2.00  where  the  effects  of  boundary- layer  and  shock-wave 
interaction  are  large.  The  explanation  for  the  independence  fr.om  a is 
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that  the  higher-order  effects  on  the  top  surface  are  compensated  for  by 
higher-order  effects  of  the  same  magnitude  on  the  lower  surface  so  that 
the  net  loading  per  unit  angle  is  very  nearly  -independent  of  angle  of 
attack. 

The  theoretical  span  loading  distribution  shown  in  figure  6(b)  can 
be  applied  to  wings  of  higher  aspect  ratio  by  simply  translating  the  tip 
solution  shown  beyond  point  2 to  the  tip  of  the  wing.  The  span  loading 
will  have  the  two-dimensional  value  of  a wing  alone  between  points  1 
and  2.  When  the  interference  from  the  body  and  wing  tip  overlaps  as  in 
figure  6(a),  the  effects  of  these  two  factors  must  first  be  separated 
before  translating  the  tip  solution. 


Effect  of  Wing-Incidence  Angle 


Comparison  is  made  between  theory  and  experiment  for  the  wing- 
incidence  case  for  data  taken  at  a Reynolds  number  of  I.5XIO®  and  Mach 
numbers  of  1.48  and  2.00  with  a = 0°  and  i^  = -1.9*^  and  It  will 

be  remembered  from  the  section  on  the  accuracy  of  data  that  there  is  no 
appreciable  error  due  to  stream  angle  for  the  wing-incidence  case. 

Pressure  distribution  in  wing -body  Junctures.-  The  linear  theory  and 
experimental  pressure  distributions  in  the  wing-body  juncture  are  compared 
in  figures  7(a)  and  7(b) . The  symbols  in  the  figures  are  flagged  to  be 
consistent  with  the  use  of  flagged  symbols  for  negative  angle-of -attack 
data.  The  figures  show  that  the  experimental  values  are  about  5 percent 
below  those  predicted  by  the  theory  for  iy  = -1.9°*  The  magnitude  of  the 
nonlinear  effects  due  to  i^^  is  predicted  at  the  leading  edge  by  shock- 
expansion  theory.  Figures  7(a)  and  7(b)  show  that  the  spread  predicted 
in  this  manner  is  in  good  accord  with  the  experimental  results.  The  pre- 
mature increase  in  the  magnitude  of  PP/ly  near  point  1 is  due  to  the 
effect  of  the  opposite  wing  panel  and  variation  of  the  local  Mach  number 
as  discussed  in  the  angle -of -attack  section.  No  significant  effect  of 
Mach  niomber  was  found  on  the  parameter  pp/iy. 

Pressure  distribution  on  top  meridian  of  body  of  wing-body  combina- 
tion.- A comparison  between  the  linear  theory  and  experiment  for  the  pres- 
sure distribution  on  the  top  meridian  of  the  body  is  made  in  figures  7(c) 
and  7(d) . These  figures  show  that  theory  and  experiment  are  in  good  accord 
for  1;;  = -1.9° • However,  nonlinear  effects  due  to  iy  cause  much  larger 
differences  between  theory  and  experiment  for  iy  = -5.7°*  This  is  consis- 
tent with  the  angle -of -attack  case  where  the  higher-order  effect  due  to  a 
was  large  for  negative  angles  of  attack. 

All  of  the  effects  observed  for  the  angle-of-attack  case  due  to  dis- 
turbances from  the  wing  are  also  shown  to  occur  for  the  wing-incidence  case 
in  figures  7(c)  and  7(d).  The  paths  of  these  disturbances  as  predicted  by 
linear  theory  are  shown  on  the  sketch  and  the  positions  at  which  the 
effects  are  expected  to  occur  are  shown  on  the  abscissa. 
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The  only  significant  effect  of  Mach  number  apparent  in  figures  7(c) 
and  7(d)  is  the  larger  boundary-layer  and  shock-wave  interaction  for 
M = 2.00  than  for  M = 1.48  near  point  1.  The  M = 2.00  experimental 
data  for  iy  = -rl.9°  dip  and  then  overshoot  at  this  point.  This  phenom- 
enon is  discussed  in  more  detail  in  the  section  of  the  report  on  Reynolds 
number  effect. 

Pressure  distribution  on  0 = 4^*^  meridian  of  body  of  wing -body  combi  - 
nation.-  Linear  theory  is  compared  with  experimental  results  for  the  pres- 
sure distribution  on  the  9 = 45°  meridian  of  the  body  of  the  combination 
in  figures  7(e)  and  7(f).  The  effects  shown  by  the  figure  are  consistent 
with  those  shown  for  the  angle -of -attack  case  and  for  the  wing-incidence 
case  on  the  top  meridian  of  the  body.  No  unusual  effects  are  observed. 

Pressure  distribution  on  wing  of  wing-body  combination.-  A comparison 
between  linear  theory  and  experiment  for  the  pressure  distribution  along 
several  spanwise  stations  is  made  in  figures  7(g)  to  7(n) . The  theoreti- 
cal values  behind  the  Mach  wave  were  obtained  directly  from  reference  7* 
These  values  were  calculated  by  essentially  the  same  method  used  in  the 
present  report,  but  without  the  aid  of  the  W2n(x,r)  functions.  The 
experimental  data  (figs.  7(k)  and  7(l))  show  that,  in  general,  pP/ly  for 
the  iy  = -1.9°  data  is  constant  and  nearly  equal  to  -2  in  front  of  the 
Mach  cone . Behind  the  Mach  cone  the  theory  generally  predicts  values 
about  5 percent  above  the  experimental  data  for  i-^/  = -1.9°  • The  higher- 
order  effects  due  to  iy  cause  larger  differences  between  linear  theory 
and  experiment  for  i^  = -5*7°*  The  figures  show  that  these  differences 
are  well  predicted  by  shock-expansion  theory.  The  effects  due  to  the 
influence  of  the  Mach  waves  are  the  same  as  those  discussed  for  the  angle- 
of -attack  case.  There  is  no  effect  of  Mach  number  evident  on  the  wing  of 
the  wing-body  combination  other  than  that  predicted  by  linear  theory. 

Span  loading  distribution.-  A comparison  between  the  theoretical  and 
experimental  results  for  span  loading  distribution  on  the  wing  and  body 
of  the  combination  is  made  in  figure  8 for  iy  = -1.9°*  The  decrease  in 
the  span  loading  due  to  the  wing  tip  was  calculated  by  the  method  of 
Busemann  (ref.  13)*  In  part  (a)  of  figure  8,  interference  from  both  the 
body  and  the  wing  tip  is  felt  between  points  1 and  2,  but  in  part  (b)  no 
interference  is  felt  between  points  1 and  2,  and  the  span  loading  is  that 
of  a two-dimensional  wing  alone. 

Figure  8 shows  that,  in  general,  the  experiment  is  5 percent  lower 
than  the  linear-theory  prediction.  Since  all  pressure  measurements  for 
the  wing-incidence  case  were  made  for  negative  values  of  iy,  the  experi- 
mental values  used  in  this  figure  were  obtained  by  doubling  the  values 
of  PP/iy  obtained  for  ly  = -1.9°  rather  than  by  considering  two  sur- 
faces as  for  the  angle -of -attack  case.  Since  this  increases  the  nonlinear 
effects  of  ly  rather  than  minimizing  them,  only  the  iy  = -1.9°  data 
(for  which  the  nonlinear  effects  are  small)  were  plotted.  However,  the 
method  of  references  6 and  7 Is  applicable  to  the  prediction  of  the  net 
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span  loading  for  larger  values  of  ly  because  the  nonlinear  effects  on 
the  upper  and  lower  surfaces  tend  to  cancel  each  other,  as  shown  for  the 
angle -of -attack  case. 


Effect  of  Reynolds  Number 


The  primary  effect  of  Reynolds  number  in  this  Investigation  was  on 
the  body.  Reynolds  n\amber  was  found  to  have  no  significant  effect  on  the 
pressure  distribution  on  the  wing  of  the  combination  for  the  range  inves- 
tigated. Figure  9 shows  the  boundary- layer  condition,  as  observed  in 
schlieren  pictures,  on  top  of  the  body  at  the  point  of  intersection  with 
the  Mach  wave  from  the  leading  edge  of  the  wing -body  Juncture  for  R = 0.6 
and  1.5x10®.  The  transition  and  separation  regions  shown  in  figure  9 
indicate  approximately  the  ranges  of  a and  iy  in  which  the  boundary 
layer  changes  from  laminar  to  turbulent  or  separated  flow  at  the  Mach  wave 
from  the  wing -body  Juncture.  In  laminar  and  turbulent  regions  the  flow 
remains  laminar  or  turbulent  across  the  Mach  wave.  Some  of  the  Reynolds 
number  effect  shown  by  figure  9 “la-y  be  due  to  changes  in  the  turbulence 
level  of  the  wind  tunnel. 

It  is  to  be  expected  that  data  obtained  for  several  angle  combina- 
tions within  any  one  of  the  regions  shown  in  figure  9 would  show  no  sig- 
nificant differences  due  to  viscous  effects,  but  that  these  data  would 
differ  from  data  in  other  regions.  For  example,  for  M = 1.48  and 
R = 0.6x10®  the  data  for  a = -2°  with  iy  = 0°  should  differ  from  the 
data  for  a = -6°  with  iy  = 0°  because  transition  occurs  at  the  shock 
wave  for  the  latter  case  but  not  for  the  former.  That  there  is  a dif- 
ference is  shown  in  figure  10  where  the  pressure  distributions  on  top  of 
the  body  for  these  two  conditions  are  compared.  In  front  of  the  shock 
wave  the  flow  is  laminar  for  both  angles  of  attack  so  that  there  is  no 
difference  in  the  two  sets  of  data.  However,  for  a =■  -6°  transition 
occurs  at  the  shock  wave  and  the  pressure  rises  as  predicted,  while  for 
a = -2°  laminar  flow  persists  behind  the  point  at  which  the  shock  is 
expected  and  the  pressure  rise  occurs  much  later  than  predicted.  In  fact, 
the  pressure  rise  does  not  occur  until  the  transition  point  shown  in  the 
figure  is  reached  and  then  it  tends  to  overshoot.  This  phenomenon  of  the 
delayed  pressure  rise  was  observed  to  occur  whenever  laminar  flow  per- 
sisted beyond  the  point  at  which  a shock  wave  from  the  wing  was  predicted 
to  exist.  When  the  disturbance  from  the  wing  is  an  expansion  wave,  the 
pressure-coefficient  curves  rise  approximately  as  predicted,  regardless 
of  the  type  of  boundary  layer.  The  conditions  for  which  this  delayed 
pressure  rise  was  observed  to  occur  are  shown  by  the  dotted  areas  in 
figure  9-  Two  other  examples  of  this  phenomenon  may  be  seen  near  points 
1 of  figures  5(d)  and  7(d)  for  a = -2°,  iy  = 0°  and  a = 0°,  iy  = -2°, 
respectively. 
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In  figure  11,  the  pressure  distributions  on  top  of  the  body  are  com- 
pared for  three  Reynolds  numbers.  It  is  shown  that  data  for  the  two  high- 
est Reynolds  numbers,  R = 1.2  and  1.5x10®,  agree  almost  identically,  while 
the  data  for  the  lowest  Reynolds  number  differ  considerably. 


Comparison  With  Theory  From  Other  Sources 


The  three  theories  for  which  numerical  results  are  available  are  com- 
pared in  figure  12 . The  theory  due  to  Ferrari  was  obtained  by  cross- 
plotting from  a figure  in  reference  8 so  that  the  curve  shown  is  only 
approximate.  The  theoretical  curve  due  to  Morikawa  is  obtained  from  tabu- 
lated results  given  in  reference  The  experimental  data  region  was 
determined  by  the  extreme  values  obtained  for  a = ±2°  for  Mach  numbers 
1.45  and  2.00.  From  this  figure  it  appears  that  either  the  theory  of 
Morikawa  or  the  theory  of  Nielsen  may  be  used  to  predict  the  pressure  dis- 
tribution on  the  wing  of  a wing-body  combination.  Ferrari 's  theory  pre- 
dicts values  that  are  somewhat  low  at  the  leading  edge  of  the  wing,  but 
it  appears  that  if  numerical  results  were  available  beyond  x/3a  = 0.7, 
they  would  lie  within  the  experimental  range.  For  a more  complete  com- 
parison of  the  theories  of  Ferrari  and  Nielsen,  see  references  1 and  2, 
Except  for  Nielsen's  theory,  no  numerical  results  for  the  pressure  distri- 
bution on  the  body  were  available  for  comparison. 


CONCLUSIONS 


An  experimental  investigation  was  ma.de  of  the  effects  of  angle  of 
attack,  a,  wing-incidence  angle,  iy,  Mach  number,  and  Reynolds  number  on 
the  pressure  distribution  on  a rectangular -wing  and  cylindrical-body  com- 
bination. The  Mach  numbers  were  1.4-8  and  2.00;  the  Reynolds  munbers  were 

0.6,  1.2,  and  I.5XIO®;  a ranged  from  +6^  to -6^  in  2^  increments;  and  ly 
ranged  from  0°  to  -5-7°  in  1.9°  increments.  On  the  basis  of  comparisons 
made  between  the  results  of  this  investigation  and  the  theory  of  NACA 
TN  2677,  the  following  conclusions  can  be  drawn: 

1.  The  theory  of  NACA  TN  2677^  in  general,  predicts  the  pressure 
coefficients  within  tlO  percent  for  values  of  a and  iy  between  +2°  and 
-2°.  For  the  wing -incidence  case,  a = 0°,  the  theory  predicts  magnitudes 
of  pressure  coefficient  above  the  experimental  values.  For  the  angle-of- 
attack  case,  ly  = 0°,  the  theory  predicts  magnitudes  of  pressure  coeffi- 
cient below  the  experimental  values. 

2.  Nonlinear  effects  due  to  angle  of  attack  and  wing -incidence  angle 
are  large.  On  the  wing  the  difference  from  linear  theory  due  to  nonlinear 
effects  of  angle  can  be  predicted  by  shock -expans ion  theory,  except  near 
the  wing -body  Juncture  for  the  angle -of -attack  case. 
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3.  Span  loading  was  shown  to  be  predicted  within  ±10  percent  for 
the  wing-incidence  case  for  all  measured  values  of  ly  on  both  the  body 
and  the  wing.  Where  comparison  was  possible,  span  loading  in  the  angle-of- 
attack  case  was  also  found  to  be  predicted  within  ±10  percent.  The  pre- 
dicted span  loadings  are  high  for  the  wing-incidence  case  and  low  for  the 
angle -of -attack  case. 

4.  For  the  angle-of -attack  case,  the  pressure  coefficients  on  the 
wing  are  experimentally  about  5 percent  higher  for  M = 2.00  than  for 
M = 1.48,  when  reduced  to  a form  that  is  theoretically  independent  of 
Mach  number.  Otherwise  Mach  n\amber  has  no  important  effect. 

5.  Viscous  effects  are  important  only  on  the  body  where  the  shock 
wave  from  the  wing  causes  large  boundary-layer  and  shock-wave  interactions 
for  some  angle  conditions . 

Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Jan.  4,  1954 
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APPENDIX  A 

DETERMINATION  OF  WING-ALONE  POTENTIAL  FOR  ANGLE -OF -ATTACK  CASE 


The  first  step  in  determining  the  potential  for  the  wing  alone  is  to 
set  up  a mathematical  model.  Since  the  exposed  wing  of  the  combination 
operates  in  the  body  upwash  field  which  effectively  twists  the  wing,  the 
wing-alone  model®  is  considered  to  be  twisted  in  the  manner  predicted  by 
Beskin  upwash  theory  for  y > a 


The  concealed  wing  may  be  extended  through  the  body  region  in  any  manner 
but,  since  equation  (Al)  gives  oy  = 2a  at  both  wing -body  Junctures,  it 
is  taken  as  a flat  plate  at  angle  of  attack  2a  (see  fig.  13). 

The  twisting  of  the  wing  is  accomplished  by  superimposing  a series  of 
flat-plate  wings  upon  a basic  flat-plate  wing  at  = 2a  (see  fig.  ll4-) . 
Each  of  the  superimposed  wings  is  at  an  incremental  angle  of  attack  and 
each  successive  wing  terminates  at  an  incremental  value  of  y greater  than 
the  previous  one.  As  these  incremental  values  of  a and  y become  infini- 
tesimally small,  the  resulting  potential  approaches  that  of  a wing  with 
the  twist  defined  by  equation  (Al) . 

For  the  purposes  of  determining  the  wing-alone  potential,  the  wing  is 
considered  to  be  composed  of  the  three  parts  shown  in  figure  13:  the  right 

exposed  half -wing,  the  left  exposed  half -wing,  and  the  wing  section  inside 
the  body.  The  perturbation  velocity  potential  is  determined  for  each  of 
these  wing  sections  and  the  results  added  together  to  obtain  the  potential 
for  the  entire  wing  alone.  Thus, 

'PW  = 'PWp  + 'PWj^  + 'PWg  (A2) 

Since  the  wing  may  be  considered  to  be  composed  of  an  infinite  niimber  of 
flat,  rectangular  wings,  the  expression 


9(aw»y)=  ^ X 


V 


cos 


- y cosh 


-1 


4^. 


+ Z COS“l 


-xy 


■+Z 


(A3) 


^oth  a and  p are  taken  as  unity.  See  Theoretical  Considerations  section 
of  text  for  method  of  generalization  to  other  values  of  these  parameters. 
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from  reference  l4  for  the  velocity  potential  of  a flat,  rectangular  wing 
will  be  used  as  the  basic  relation  for  the  calculations.  Equation  (A3) 
gives  the  velocity  potential  at  any  point  (x,y,z)  due  to  a flat,  rectangu- 
lar wing  at  angle  of  attack  terminating  at  y = 0,  and  extending  to 

•0  along  the  positive  y axis. 

Since  the  twisted  wing  was  shown  to  be  equivalent  to  a basic  flat- 
plate  wing  at  angle  of  attack  2a  plus  an  infinite  number  of  modifying 
flat-plate  wings,  the  potential  of  the  right  half -wing  may  be  written  as 

= (p(2a,y-l) +£iCp(Axi,y-T1i)  (a4) 

R 

The  first  term  on  the  right  in  equation  (a4-)  is  the  potential  due  to  the 
basic  exposed  half -wing  which  terminates  at  the  wing -body  juncture,  y = 1, 
and  is  at  angle  of  attack  2a.  The  second  term  is  the  potential  of  the  i 
modifying  wings  each  at  angle  of  attack  Aa^  and  terminating  at  y = 
where  1 < t)^<  «.  Since  equation  (A3)  is  homogeneous  with  respect  to  a, 
equation  (a4)  may  be  written 

<p„  = (p(2a,y-l)+ 2iCp(l,y-Tii)Axi  (A5) 

From  equation  (Al) 

= - ft  (a6) 

Therefore, 

= <p(2a,y-l)  - 2aJ  (A?) 

where  the  limits  of  integration  are  determined  by  the  range  of  y on  the 
wing  included  in  the  fore  Mach  cone  originating  from  the.  point  for  which 

is  being  determined.  From  figure  15  it  is  apparent  from  the  equa- 
^E 

tion  of  the  fore  Mach  cone  emanating  from  the  arbitrary  point  Pi  for 
p = 1 that 

= zi^  + (Tiu-yi)^ 

Therefore,  the  upper  limit  of  integration  is 

Tlu  = 7i  + ^ Xi^-  zp 

The  lower  limit  of  integration  is  at  the  wing-body  juncture,  y 
equations  (A3);  (A7),  and  (a8) 


(a8) 

= 1 . From 
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- X cos 


-1  -(y-1) 


(y-l)  cosh 


-1 


*/  (y-i)^+z^ 


z cos 


-1 


- x(y-l) 


2VoC£. 


v/(y-l)^+z2  a/x^-z^  J 

(y-T])  cosh"^ 


^ (y-n) 


+ z cos 


-1 


y+-\/  x2_z2 


-x(y-T)) 


X cos 


-(y-nl 


2+z2 


(A9) 


Carrying  out  the  integration  and  combining  terms  gives 

y 


Va 

It 


- X cos“^ — +1  ^ 


( ^ - y) 


^/(y-1) 


2+z2 


y — oo»h-i  + 

(y^+z^)  a/x^ -y^  -z^  y x^-z^ 


z + cos*^ 


y2+z2y/  y(y-l)2tz2  V^X2-Z2 


]} 


(AlO) 


Equation  (AlO)  gives  the  potential  due  to  the  exposed  right  half -wing.  To 
this  must  be  added  the  potentials  due  to  the  other  two  wing  sections.  The 
potential  due  to  the  section  of  the  wing  in  the  body  region  is  simply  the 
difference  between  the  potentials  of  two  flat  wings  at  Oy  = 2a.  One  of 
these  wings  terminates  at  the  wing-body  Juncture  at  y = -1,  figure  13^ 
and  extends  (through  the  body)  indefinitely  in  the  positive  direction. 

The  other  wing  terminates  at  the  other  wing -body  Juncture,  y = +1,  and 
also  extends  indefinitely  in  th6  positive  y direction.  The  difference 
between  the  potentials  of  these  two  wings  is  the  potential  of  the  wing 
section  in  the  body 


=»  <p(2a,y+l)  - 9(2a,y-l)  (All) 

B 

The  expression  for  9(o^j,y)  is  given  by  equation  (A3) . Since  the  model 
is  symmetrical  about  the  vertical  y = 0 plane,  the  potential  for  the 
other  half -wing  is  simply  obtained  by  replacing  y by  -y  in  equation 
(AlO) . 


(x^-y»z)  (A12) 

Combining  equations  (AlO),  (All),  and  (A12)  gives  the  potential  due  to  the 
entire  wing  alone. 
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(AI3) 


Investigation  of  <Pw  given  by  equation  (A13)  reveals  that  there 
are  three  regions  on  the  body  in  which  the  real  part  of  this  expression 
differs.  A fourth  region,  region  IV,  is  entirely  on  the  wing  and  is, 
therefore,  not  necessary  for  determining  the  normal  velocity  distribution 
on  the  body.  These  regions  are  determined  by  three  characteristic  Mach 
cones.  One  of  the  Mach  cones  originates  on  the  body  axis  at  x = 0,  and 
the  other  two  originate  at  the  leading  edges  of  the  two  wing -body  Junc- 
tures (see  fig.  16).  The  expressions  for  the  real  part  of  in  the 

three  regions  on  the  body  are: 

Region  I: 
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Region  II: 

CD  - Va 

'Pw  - - V 


-X  coB-i- + y-2')  cosh-i 
^/x2-z2  \y^+Z^  / 


V(y-l)^+z" 


/ z \ _,  - x(y-l) 

( - Z I COS  1 - ' , 

y2+z2  J y(y_l)2^.22  Vx2-z2 


J2L 


cosh 


(y2+z^)  -/x2-y2-z2 


_i  x^-z^-y(y-l).  2 (z-x)/ 


Region  III; 
Va 

^Pw  = - - 


-X  cos 


-1  ^ y + y-2^  cosh“i 

Vy^+^2  y 


^(y-i) 
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z cos 
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APPENDIX  B 

DETERMINATION  OF  VELOCITY  AMPLITUDE  FUNCTIONS 


The  method  of  references  6 and  7 requires  a Fourier  series  expansion 
of  the  expression  for  the  normal  velocity  on  the  body  surface.  The  coef- 
ficients, fan(^)^  cosine  terms  in  this  series  are  called  the 

velocity  amplitude  functions.  The  interference  pressure  coefficients  are 
found  directly  from  these  functions,  as  will  be  shown  in  Appendix  C. 


The  velocity  normal  to  the  body  surface  due  to  the  wing  alone  is 
easily  found  from  the  wing-alone  velocity  potential  by  differentiation. 
The  amplitudes  of  the  Fourier  series  expansion  of  the  normal  velocity 
expression  are  obtained  by  the  usual  equation, 


sin" 


dC 


^2n 


sin 


-1 


cos(2n0)d0 


^2n(*)  functions  were  obtained  in  closed  form  for  the  wing-incidence 
case,  but  the  angle -of -attack  case  required  the  evaluation  of  untabulated 
elliptic  integrals  so  that  numerical  methods  had  to  be  used  to  obtain 
f2n(x).  The  resulting  f2n(x)  functions  are  shown  in  figure  I7.  This 


figure  shows  that  the  curves  of 

Va 

between  x =■  0 and  x = 1 but  for  x>l, 
is  nearly  constant. 


against 

(x) 


‘■zn 


Va 


X oscillate  n times 
ceases  to  oscillate  and 
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APPENDIX  C 

PRESSURE  COEFFICIENTS  OF  THE  WING-BODY  COMBINATION 
FOR  THE  ANGLE-OF-ATTACK  CASE 


Determination  of  Interference  Pressure  Coefficients 


The  interference  pressure  coefficients  ar6  obtained  directly  from 
the  velocity  amplitude  functions,  f2n(x)  f using  the  universal  functions 
W2n(x,r),  developed  in  references  6 and  7*  In  references  6 and  7 it  is 
shown  that  the  interference  pressure  coefficients  for  r = 1 for  each  of 
the  fan(^)  functions  are 


^*2n  2n0 


-/ 


f2n(x)  -/  f2nU)W2n(x-IA)dS 


(Cl) 


The  integration  operation  requires  graphical  or  numerical  methods.  The 
total  interference  pressure  coefficient  is  the  svim  of  the  n components, 
P2n.  To  get  the  exact  linear-theory  solution  it  would  be  necessary  to 
compute  an  infinite  niimber  of  terms,  but  in  most  cases  using  only  four 
terms  gives  a close  approximation  to  linear  theory,  as  shown  in  reference 
7 and  verified  by  figures  18,  19,  and  20.  A relation  similar  to  that  of 
equation  (Cl)  gives  values  of  ^2n  for  r different  from  1. 


Pressure  Distribution  in  Juncture  of  Wing-Body  Combination 


The  pressure  distribution  of  the  combination  is  obtained  by  adding 
the  interference  pressure  coefficients  to  the  pressure  coefficients  of 
the  wing  along.  The  results,  using  four  and  six  Ppn  components,  are 
shown  in  figure  I8.  This  figure  shows  that  four  components  give  a close 
approximation  to  the  linear-theory  value  for  x/pa  > 1.  At  x/pa  = 0, 
the  wing  leading  edge,  linear  theory  using  Beskin  upwash  gives  exactly 
pp/a  = -4.0.  For  the  region  x/pa  < 1 the  higher  harmonics  have  their 
greatest  importance  and  many  components  would  be  necessary  to  get  good 
accuracy.  However,  satisfactory  accuracy  can  be  obtained  by  fairing  a 
curve  through  this  region  since  both  end  points  are  known. 

One  point  of  Interest  in  figure  18  is  the  increase  in  the  magnitude 
of  pP/a  near  point  1.  This  is  due  to  the  Influence  of  the  opposite  half- 
wing  reaching  the  wing -body  Juncture  at  this  point  as  shown  in  the  sketch. 


Pressure  Distribution  on  the  Body  of  the  Wing-Body  Combination 


The  pressure  distribution  on  the  body  is  also  obtained  by  adding 
the  interference  pressure  coefficients  to  the  pressure  coefficients  due 
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to  the  wing  alone.  The  interference  pressure  distribution  for  any  value 
of  6 differs  from  that  in  the  wing -body  Juncture,  6 = 0,  only  by  the 
cos  2nd  factor  in  equation  (Cl) . For  example,  in  the  juncture,  cos  2n& 
is  always  +1,  on  top  of  the  body,  6 = jt/2,  cos  2n0  alternates  between 
+1  and  -1  as  n increases,  and  on  the  6 = n/k  meridian  cos  2n0  has 
values  of  0,  +1,  and  -1  so  that  when  n is  odd,  ^2n  = 0.  The  pressure 
distributions  on  the  top  meridian  of  the  body  and  on  the  0 = 45°  merid- 
ian of  the  body  are  shown  in  figures  19  and  20,  respectively. 

Several  interesting  effects  are  esAibited  by  figures  I9  and  20.  The 
step  in  the  wing-alone  pressure  at  x/pa  = 1 in  figure  19  is  effectively 
canceled  by  the  interference  pressure  from  x/pa  = 1 to  x/pa  = n/2  and 
for  x/pa  > jt/2  the  pressure  increases  rapidly.  The  effect  of  the  inter- 
ference pressure  in  canceling  the  effect  of  the  wing  alone  on  the  top  of 
the  body  from  x/pa  = 1 to  x/pa  = n/2  is  to  be  expected  since  the  wing 
of  the  combination  can  have  no  effect  on  the  body  in  front  of  the  Mach 
helix  (point  1 of  sketch)  originating  at  the  leading  edge  of  the  wing-body 
Jimcture.  If  an  infinite  number  of  components  had  been  computed,  the  com- 
bination pressure  coefficients  would  be  identically  zero  from  x/pa  = 0 
to  x/pa  = n/2.  The  same  effects  are  exhibited  by  figure  20  except  that 
the  wing-alone  step  occurs  at  x/pa  = ,j2/2  and  the  Mach  helix  intersects 
the  meridian  at  x/pa  = jt/4,  point  1.  The  Mach  helix  from  the  opposite 
wing  panel  intersects  the  meridian  at  point  2 causing  an  additional  pres- 
sure rise.  Since  the  region  in  which  pp/a  = 0 is  known  and  since  the 
e^cact  linear  theory  is  well  approximated  by  four  components  for  large 
values  of  x/pa,  theoretical  curves  of  good  accuracy  can  be  faired  from 
figures  19  and  20.  The  area  under  the  high  peaks  in  the  curves  near 
x/pa  = jt/4  would  become  infinitesimal  if  an  infinite  number  of  inter- 
ference pressure  components  were  taken. 


Pressure  Distributions  for  the  Wing  of  Wing-Body  Combinations 


In  region  IV  of  figure  16,  the  calculation  of  pressure  coefficients 
is  just  a wing-alone  problem.  The  pressure  coefficients  in  this  region 
can  therefore  be  obtained  directly  from  the  wing-alone  potential  as  given 
by  equation  (AlO).  The  result  is 


P = 


2a 


1 + 


y 

(y2.x2)3/2 


(C2) 


In  the  region  behind  the  Mach  cone  the  pressure  coefficients  can  be 
obtained  directly  from  the  functions,  as  was  done  on  the  body. 

The  span  loading  distribution  is  determined  by  graphical  integration 
of  faired  pressure-distribution  curves  of  figures  18,  19,  and  20.  The 
pressure  results  of  these  figures  are  for  c/pa  < 4 and  for  pA  > 2 . Span 
loadings  for  any  combination  of  c/pa  and  PA  in  these  ranges  can  be 
obtained  from  these  pressure  distributions.  The  theoretical  and  experi- 
mental span  loading  distributions  on  the  body  were  evaluated  back  to 
the  wing  trailing  edge. 
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For  purposes  of  specifying  the  span  loading,  the  following  equation 
is  used: 


> (C3) 


The  span  loading  is  the  quantity  inside  the  bracket. 
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TABLE  I.-  ORIFICE  LOCATIONS  ON  WING  AND  BODY  OF  WING-BODY  COMBINATION 
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Two  independent  measurements 


juncture  for  a = 0°,  = -1. 9°,  M - 1. 48,  and  R - 1. 5x  10®. 
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(g)  M=  1.48,  y/a  = 1.25 


4 


-3 


/ 


0 

.5  10  L5  2.0  25  3.0  3.5 


x//3a 


(h)  M=  2.00,  y/a  = ,1.25 
Figure  7.-  Continued. 


NACA  TN  3128 

Laminar 

WWVs!  Transition 

[III III  Turbulent 

V/////  Separation 

iliili  Delayed  pressure  rise  region 

(a)  M = 1.48,  R=  0.6x10® 


(b)  M = 1.48,  R = 1.5x10 


Wmg  -incidence  angla,  deg 

(c)M  = 2.00,  R = 0.6X10®  (d)M=2.00,  R =1.5x10® 

Figure  9.-  Boundary- layer  condition  on  top  meridian  of  body  at  point  of 
intersection  with  Mach  helix  from  wing. 
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Untar  theory- 
TN  2677 


Transition  Transi  tion 

point 

\fora:-2*\ 

L6  24  32  4.0  4.8 

x/fia 

Figure  10.-  Effect  of  transition  position  on  pressure  distribution  on  top 
meridian  of  body  with  M = 1. 48,  R = 0. 6x10®,  and  i'^  = 0°. 


for  a=-6 


Linear  theory- 
TN  2677 


I / er  RU.5xtO\ 

Ef  1.2x10^ 

I ^ 0.6x10^ 

0 .8  t:6  2.4  32  4.0  4.8 

X/0<S 

Figure  11.-  Effect  of  Reynolds  number  oh  pressure  distribution  on  top 
meridian  of  body  with  M •=  1.48,  .a  = -6°,  »td  i^  = 0°. 
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I 


Figure  13.-  Shape  of  wing  alone  with  effective  twist  produced  by  body 
upwash  field  for  angle-of -attack  case . 


Basic  wing 


Twisted  wing 


Figure  14.  - Formation  of  twisted  wing  by  superposition  of  Infinite 
number  of  fiat  plates. 


Figure  18.-  Theoretical  pressure  distribution  in  wing-body  juncture  of  combination 
approximated  by  use  of  four  and  six  Fourier  components. 
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NACA-Langlejr  - 4-19-54  - 1000 


Figure  20.  - Theoretical  pressure  distribution  on  9 = 45°  meridian  of  body  of  combination 
approximated  by  use  of  four  and  six  Fourier  components. 
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